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We examine whether or not the surface tension acts to increase the nucleon density in the nuclear 
interior within a compressible liquid-drop model. We find that it depends on the density dependence 
of the surface tension, which may in turn be deduced from the neutron skin thickness of stable nuclei. 



o 
o 



in 
o 

o 

o 

> 

I 



X 



PACS numbers: 21.10.Gv, 21.65.+f 



The saturation property of bulk nuclear matter is usu- 
ally deduced from empirical data on the masses and radii 
of stable nuclei 1]. If nuclear matter is incompressible, 
the bulk saturation density is equal to the density in the 
nuclear interior. However, nuclear matter in nuclei is 
more or less compressible. The question of how the finite 
compressibility of nuclear matter makes the density in 
the nuclear interior deviate from the saturation density 
has yet to be clarified. In this paper, by using a com- 
pressible liquid-drop model, we show that a key feature 
in addressing this question is the density dependence of 
the surface tension. We then derive its relation with the 
neutron skin thickness. 

The dependence of the surface energy on the density 
of the nucleon liquid naturally occurs. In the Fermi-gas 
model, the surface energy arises from a reduction of the 
available density of states of nucleons due to the presence 
of the surface, beyond which the nucleon flux vanishes. 
The total kinetic energy thus increases by an amount 
proportional to the surface area, and this increase de- 
pends on the nucleon density through the nucleon Fermi 
momentum. In addition to this kinetic contribution, the 
interaction between nucleons contributes to the surface 
energy more importantly. This is because a nucleon in 
the surface region does not perceive the same amount of 
attraction from the surrounding nucleons as that which it 
would if it were in the deeper region. This contribution is 
determined by the property of the nuclear medium, which 
is characterized by the nucleon density. We remark that 
for the same reasons, the surface energy depends also on 
the neutron excess of the liquid; this dependence has to 
be invariant under exchange between neutrons and pro- 
tons. 

As discussed by Yamada 0, the dependence of the 
surface energy on the inner liquid density controls the 
density deviation from the bulk saturation density. Usu- 
ally, this dependence is not considered explicitly. This is 
because the surface tension is normally calculated for a 
planar interface between the saturated nucleon liquid and 
the vacuum, which are in mechanical equilibrium. How- 
ever, this is not equivalent to mechanical equilibrium in a 
real nucleus, which is generally associated with additional 
pressures arising from the size- and density-dependence 



of the surface energy. These pressures in turn affect the 
equilibrium density of the compressible nucleon liquid. 

In order to see this effect, we utilize a compressible 
liquid-drop model, which gives rise to a semi-empirical 
mass formula in a way dependent on the density and neu- 
tron excess in the nuclear interior, nin and Jin- Generally, 
a liquid-drop model is advantageous to the description of 
various macroscopic properties of nuclei. During the past 
decade, it has been used in describing, e.g., neutron skin 
13 , nuclear fission , deformation of rapidly rotating nu- 
clei 0, synthesis of superheavy nuclei Q, and nuclei in 
neutron star crusts 0. 

Throughout this paper we consider nearly symmetric 
nuclei, for which we can set i?„ ~ Rp. We assume that 
the distribution of i nucleons (i = n,p) is spherically 
symmetric, uniform at a number density n^, and squared 
off at a radius Ri . For a nucleus of mass number A and 
charge number Z (neutron number = ^ — Z), we thus 



obtain 



and 



Sin = [nn - np)/ (n„ -|- rip) 



{N -Z)/A. 



(1) 



We then write the binding energy Eb of the nucleus as 

— Eb = Evol + -Esurf + -Bcoul ■ (2) 

Here 

£^voi = ^w(nin,(5i„), (3) 

with the bulk energy per nucleon w, is the volume energy, 

Esurf = 4:TTa{nin,Sin)Rl, (4) 

with the density-dependent surface tension a, is the sur- 
face energy, and 



5i?„ 



(5) 



is the Coulomb energy. For w and a, we adopt a form 
expanded with respect to the density and neutron excess 
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around = uq and 5-^ = 0: 



Here 



w{nin,Sin) = Wo + -^^(riin - UoY 



(6) 



where no and wq are the saturation density and energy 
of symmetric nuclear matter, Kq is the incompressibihty 
of symmetric nuclear matter, is the symmetry energy 
coefRcient, and L is the density symmetry coefficient, and 



cr{nin, Sin) = CTq 



1 - Csym^in + X 



riQ 



no 



(7) 



where ctq = a(no,0), Cgym is the surface symmetry en- 
ergy coefficient, and x = ino/cro)da/dnin\ni^=no,5i,-,=Q- 
Eq. ^ we have ignored the energy contribution of the 
neutron skin thickness Rn — Rp, which will be considered 
later, and curvature corrections. We have also ignored 
pairing and shell corrections since we will confine our- 
selves to macroscopic properties of the nuclear ground 
state. We remark that in equilibrium riin is related to 
Sin, as we shall see just below. 

Some of the coefficients characterizing the bulk energy 
© can be deduced from empirical data for the masses 
and root-mean-square charge radii of stable nuclei. The 
saturation density no, the saturation energy wo, and the 
symmetry energy coefficient 5*0 typically take on a value 
ranging 0.14-0.17 fm-^, -16 ± 1 MeV, and 25-40 MeV, 
whereas the incompressibihty Kq and the density symme- 
try coefficient L, which control the density dependence 
of bulk nuclear matter, are not well constrained. Using 
a simplified version of the Thomas-Fermi model we 
found that various sets of the values of Kq and L ranging 
180-360 MeV and 0-200 MeV reasonably reproduce the 
empirical masses and radii and that future systematic 
measurements of the matter radii of unstable neutron- 
rich nuclei would give a good constraint on the value of 
L. 

We turn to the coefficients in the surface tension |(7|). 
The primary coefficient cto and the surface symmetry co- 
efficient Csym can be estimated from the empirical mass 
data as CTo — 1 MeV fm~^ and Cgym = 1.5-2.5. The pa- 
rameter X characterizing the density dependence of the 
surface tension is poorly known and hence the quantity 
of interest in this work. Myers and Swiatecki Q simply 
set X = 0, while the Fermi-gas model predicts x — 4/3. 
We will see that a precise determination of no from the 
values of Uin deduced, e.g., from electron- nucleus elastic 
scattering data requires reliable information about x- 

Let us now estimate the equilibrium value of nin from 
pressure equilibrium and compare it to the bulk satu- 
ration density n^ at fixed Sm- Within the present com- 
pressible liquid-drop model, the pressure equilibrium con- 
dition can be obtained from optimization of the binding 
energy ((2Jl with respect to the size under fixed A and Z 
as 



= R 



vol 



Pi 



Coul- 



(8) 
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(9) 



is the volume pressure, 



2cto 

Rn 



1 - oX - Csyn^Sfn + X 



njn - no 
no 



is the surface pressure, and 



Po 



'oul 



3ZV 
207ri?4 ' 



(10) 
(11) 



is the Coulomb pressure. The bulk pressure vanishes at 
the saturation density, 



no 



3ino 2 
Ko 



(12) 



which generally decreases with increasing Sm as has al- 
ready been discussed in Refs. 0,0]. The Coulomb pres- 
sure acts to increase the nuclear size, whereas the surface 
pressure tends to enlarge or reduce the nuclear size ac- 
cording to whether x is larger or smaller than ~ 2/3. 

The relation (|SJ), if the Coulomb pressure Pcoui is ig- 
nored, can be reduced to Laplace's formula. This can be 
done by transforming Eq. Q into 



vol 



3croX 2ct( 



Rv 



Rv 



(13) 



Here the left side, arising from the energy derivative with 
respect to nin, corresponds to the pressure of the nucleon 
liquid, while the right side arises from the energy deriva- 
tive with respect to Rp. 

Deviation of the equilibrium value of n^n from the bulk 
saturation density n^ at fixed Sin can be estimated from 
condition © as 



0.016 



230 MeV\ 



Ko 



1 MeV fm" 



5 fni 

Rn 



407ri?3cro / 



(14) 



In this estimate we have used Psmf — ~2cro(l — 3x/2)/-Rp. 
The ratio of the Coulomb pressure to the surface pres- 
sure, 3Z2e2/407ri?3(To, is typically 0.2-0.6. We find from 
this estimate that for x—^ and 4/3, the surface pressure 
can induce about 10 % change in nin in different direc- 
tions. We thus see the role played by x in determining 
ng from empirical information about njn. 

We now proceed to show that the neutron skin thick- 
ness, which has been neglected so far, is a quantity that 
may be useful for deduction of the value of x- In do- 
ing so, as considered by Pethick and Ravenhall (lol| , it is 
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convenient to describe the nuclear surface in a thermo- 
dynamically consistent manner. In this description, the 
nuclear surface is in thermodynamic equilibrium with the 
bulk system composed of A nucleons, and the neutron 
skin arises from adsorption of Ns neutrons onto the nu- 
clear surface. The interior region composed oi A — Ns 
nucleons acts as a reservoir of neutron chemical poten- 
tial and neutrons can go back and forth between 
the skin and interior regions. Consequently, the relevant 
thermodynamic quantity is the thermodynamic poten- 
tial, fl ~ ilvoi + f^surf + ^^Couh divided in a similar way 
to the binding energy In equilibrium, flguri — crA, 
where A is the surface area. A small quasistatic change 
in the neutron excess in the interior region with A and N 
fixed gives rise to a change in the thermodynamic poten- 
tial of the surface, Ar^smf, and a change in the neutron 
chemical potential, A/i„, which are related as 



(15) 



This relation indicates that the neutron skin can be de- 
scribed in terms of the bulk and surface properties. We 
remark that at N = Z, a balance between the induced 
changes Afivoi and Aficoui in the volume and Coulomb 
energies allows the neutron excess in the interior region 
to deviate from zero and hence a proton skin to occur, 
as we shall see. 

It is straightforward to combine the above thermody- 
namic description of the nuclear surface with the com- 
pressible liquid-drop model adopted here. In this model, 
Sim -^1 a-iid /i„ read 



N -Ns- Z 

A-Ns 
N - Z 3{Rn - Rp) 



2R„ 



Ns ~ AirRlnniRn - Rp), ^ ~ 47ri?^ 



and 



Wo + 5o<5i„(2 - Si,,) + 0{Sl). 



(16) 



(17) 



(18) 



Here we have calculated ^„ at the bulk saturation density 
(|12|l : the dependence of /i„ on the parameters L and Kq 
characterizing the density dependence of the bulk energy 
does not appear up to second order in Si^,. 

We can now obtain the expression for the neutron skin 
thickness in the absence of Coulomb energy. In this case 
the system is symmetric under exchange between neu- 
trons and protons. Substitution of Eqs. (Tj), ((T7|) . 
and lfTS|l into Eq. (fTK|l leads to 



Rn — Rp — C(5 I 1 



3C 
2i?B 



+ OiS^), (19) 



where 



c = ^(c +^ 
^- Som r^^"^+ Ko 



(20) 



and S = {N — Z)/A. The parameter C originates mainly 
from a change in the surface tension due to the small qua- 
sistatic change in Sin- This change is characterized not 
only by the surface symmetry energy coefficient Cgym, 
but also by the parameter x through the dependence of 
the saturation density Us given by Eq. (|12|l on Sin- The 
term 3Lx/Kq in the right side of Eq. (|20|) is associated 
with the finite compressibility and thus vanishes in the 
incompressible limit in which Eq. H19|l reduces to the re- 
sult for Rn — Rp obtained by Pethick and Ravenhall ■ 
Note that this term does not exist in the result of Myers 
and Swiatecki ,11] who presumed x = in a compressible 
liquid-drop picture, although it can be comparable with 
Csym- The fact that 3Lx/Ko is poorly known suggests 
that one could not deduce the EOS parameters L and 
Kq from experimental data for the neutron skin thickness 
without knowing x- Consequently, it turns out that pre- 
vious investigations that attempted to relate the neutron 
skin thickness with the EOS of nuclear matter 0, 0, [3| 
do not take full account of uncertainties in the parame- 
ters X, L, and Kq. We can see from Eq. H19(l that the 
neutron skin vanishes at N ~ Z, as it should in the ab- 
sence of Coulomb energy. 

Coulomb effects ignored in Eq. (|19|l induce a thin pro- 
ton skin at N = Z through a deviation of Sin from zero 
and a polarization of the nuclear interior, as discussed 
by Myers and Swiatecki 11]. First, in order to calculate 
the deviation of Sin at N — Z, one has only to consider 
the balance between Afivoi and Aficoui- In this case an 
increment in the proton radius Rp tends to reduce the 
Coulomb energy, while it leads to a cost of the symmetry 
energy in the nuclear interior. We may thus obtain 



Ze^ 



20i?pS'o 



at N = Z, 



(21) 



which typically amounts to a small value of ^ 0.02. For 
the neutron skin thickness H19|l . this deviation effectively 
replaces S hy S — Ze^ /20RpSo In the liquid-drop 

picture, the (3 stability condition that the neutron and 
proton chemical potentials are almost equal is given by 



Sin 



3Ze" 
10i?p'5'o 



(22) 



This Sin is much larger than Ze^ /20RpSQ, which implies 
that generally stable nuclei, lying around the (3 stability 
line, have a neutron skin. Second, the polarization of 
the nuclear interior tends to reduce a difference between 
the root-mean-square radii of the neutron and proton 
density distributions, Ar„p, by redistributing nucleons 
in such a way as to deplete protons in the central region. 
This reduction amounts to y^Ze'^/70So HIl, which 
typically takes on a small value of 0.01-0.04 fm. 

By incorporating these Coulomb effects into Eq. H19|l . 
we finally obtain 



Ar 



np 



c s- 



Ze^ 



20RpSQ 
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Ze^ 



70So 
(23) 
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FIG. 1: Difference in the root-mean-square radius between 
neutrons and protons for six stable nuclei of ^ > 50. The 
squares and crosses denote the results calculated from Eq. 
I)2i-i|l for X ~ and 4/3; the other parameters are set to be 
no = 0.16 fm-^ So = 30 MeV, Kg = 230 MeV, L = 100 
MeV, cro = 1 MeV fm-^ Csym = 1.8, and Rp = 1.2A^^^ fm. 
The empirical data (dots) are taken from Refs. [l^ Il3|. 



where a factor y 3/5 arises from a difference between the 
root-mean-square and the half-density radius for the rect- 
angular distribution. Here we have ignored corrections 
due to the difference in the surface diffuseness between 
protons and neutrons. These corrections are lower order 
in A than the terms in Eq. H23|l Q and hence expected 
to be negligibly small for heavy stable nuclei of interest 
here. 

We now ask how one can obtain information about x 



from empirically deduced values of Ar„p for stable nuclei 
such as Ni and Sn isotopes. Such values can be deduced, 
e.g., from measurements of proton- and electron-nucleus 
elastic differential cross sections Ts^. It is instructive to 
compare the deduced values with the prediction by the 
present liquid-drop model, which is given by Eq. H23|l . 
Figure 1 exhibits the deduced and predicted values as 
a function of 5. When x = 0, the predicted values by 
setting the other parameters at typical values are appre- 
ciably smaller than the deduced ones. This suggests that 
X is likely to be positive. However, the magnitude of x 
remains to be clarified since x is coupled with the un- 
certain parameters L and Kq in Eq. H23() . We remark 
that a staggering of the deduced values is far larger than 
that of the predicted values. This may be partly because 
the former values were deduced by various groups using 
different models for proton elastic scattering and nucleon 
density distribution, and partly because pairing and shell 
effects are ignored in the present prediction. 

In summary we have found from a compressible liquid- 
drop model that whether or not the nucleon density in the 
nuclear interior is larger than the bulk saturation density 
depends on the density dependence of the surface tension, 
which in turn controls the neutron excess dependence of 
the neutron skin thickness. In order to deduce the density 
dependence of the surface tension and the bulk saturation 
density from the neutron skin thickness and the interior 
density, it would be useful to systematically analyze dif- 
ferential cross sections measured for proton and electron 
elastic scattering off stable nuclei. In such analysis of pro- 
ton elastic scattering data obtained for incident energies 
above 500 MeV, one could relate the angle of diffraction 
maxima measured in the small momentum transfer re- 
gion to the root-mean-square matter radius by using the 
Glauber theory in the optical limit approximation |l5l |. 
The research in this direction is under way. 
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